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Abstract 

The spin | XXZ chain in a random magnetic field pointing in the Z direction is nu- 
merically studied using the Density Matrix Renormalization Group (DMRG) method. The 
phase diagram as a function of the anisotropy of the XXZ Hamiltonian and the strength 
of the random field is analyzed by computing the spin correlations and the superfTuid den- 
sity. To obtain the superfTuid density we consider a superblock configuration representing 
a closed system with an arbitrary twist at the boundary. This allows us to estimate the 
size of the critical region where the quasi-long-range order persists, that is, where the spin 
correlation length is infinite and the superfluid density is non-zero. 
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1 Introduction 



The study of quantum phase transitions represents a wide area of research in condensed 
matter physics . These transitions occur at T = when a parameter of the hamiltonian is 
varied across a critical value and are driven by quantum fluctuations. Among them, those 
that are due to the effect of disorder are more difficult to analyze both analytically and 
numerically and have been the focus of much attention recently. 

In this paper we consider the spin 1/2 XXZ chain in the presence of a non-symmetry 
breaking random magnetic field. This system is particularly interesting since a renormal- 
ization group analysis [§, ||, |I| suggests the existence of a critical region in the space of 
parameters where the disorder is not relevant. This idea is also supported by numerical 
calculations (5|, |6| . In || the system was studied using an exact diagonalization method up 
to 16 sites. Although this size is enough to study the system far from the critical region, 
near the transition, as it is explained there, the calculations were less reliable due to the 
small size of the chains considered. Larger chains can be handled using a density matrix 
renormalization group (DMRG) algorithm [fjj at the expense of truncating the Hilbert space 
of states. This method was used in Q to obtain further information about the critical region 
by computing the phase sensitivity, which is the difference in energy between periodic and 
antiperiodic boundary conditions. 

In this paper we use the DMRG method to compute the correlation length and the 
superfluid density which give direct information about the existence of a quasi-long-range 
order region since that region is defined as the one where the correlation length diverges and 
the superfluid density is non-zero. In particular, to compute the superfluid density we used 
a superblock configuration that allowed us to introduce an arbitrary twist at the boundary. 

Since we perform DMRG calculations, our results contain two approximations. The 
first is that we consider a finite chain and the second that we truncate the Hilbert space 
(keeping a finite number of density matrix eigenstates). These approximations are related 
since with a truncated Hilbert space the correlation length is always finite and, in fact, the 
value corresponding to an infinite chain (what we are going to call the 'truncated system 
correlation length') becomes smaller as we reduce the Hilbert space. To minimize the finite 
size effects in the calculation of the spin correlations, it is necessary to consider chains 
larger than the truncated system correlation length. However, in the critical region this 
correlation length increases rapidly as we keep more DMRG states, therefore, reducing the 
truncation error by keeping a large number of states soon becomes prohibitive since we 
need to consider very long chains []. Thus, most of our calculations of the spin correlations 
are done keeping up to 32 states and give only an estimation to the shape of the critical 
region. As we will discuss later, the main source of error in this estimation is the truncation 
of the Hilbert space. To further support these qualitative results, for the special case of 
anisotropy A = —0.75 and selected values of the random magnetic field, we keep up to 
96 states and extrapolate the computed correlation length to infinite number of sites and 
infinite number of DMRG states. The extrapolations allow us to conclude the existence of 
a quasi-long-range order region and put bounds on the maximum critical field. 

In the calculations of the superfluid density we also keep up to 32 states in most of the 
cases and estimate the shape of the critical region using a criterion in which, again, the main 
approximation is the truncation of the Hilbert space. As we did for the correlations, for 

1 For a more quantitative statement see Fig. hq 
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the special case of A = —0.75 and selected values of the random magnetic field we keep up 
to 128 states and extrapolate the superfluid density to infinite length and infinite number 
of DMRG states. We find an agreement between the correlation length and the superfluid 
density results. 

The remainder of this paper is organized as follows. In section [| we summarize the 
known results for this system and in section we describe the numerical method used. We 
present and discuss the results in sections |4|and |5| and finally, in section ^ we give our 
conclusions. 

2 The spin 1/2 XXZ chain in a random z magnetic field 

In this section we summarize previous results obtained for this system by C. Doty and D. 
Fisher £§. 

In the absence of disorder the spin 1/2 XXZ chain is defined by the Hamiltonian 

H = J2( S ? S i+i + S i S i+l + ^!Sf +1 ). (1) 

i=l 

When the anisotropy parameter, A, satisfies —1 < A < 1, the system is in a non-magnetic 
gapless phase where the spin-spin correlation functions in the ground state decay as a power 
law (quasi-long-range order (quasi-LRO) phase) [§, |9|, |l0fl . On the other hand, the spectrum 
has a gap for A > 1 (A < — 1) and possesses a long-range Ising like antiferromagnetic 
(ferromagnetic) order. 

In this work we study the effects produced by the addition of the following term: 

2fi = X>*Sr, (2) 
i=l 

where /if is a random magnetic field with an independent and uniform distribution in 
each site ((hf) = 0, (hfhj) = 5ijh,Q/3). Note that the addition of this term does not break 
the symmetry in the xy plane. 

The main conclusion of 0] is that for —1 < A < —1/2 the quasi- long-range order that 
the system possesses at h\ = persists for finite values of the random magnetic field. 

To reach this conclusion, the long-distance behavior of the system was analyzed by 
studying the renormalization group flow. In the region |A| < 1 one can perform a Jordan- 
Wigner transformation to spinless fermionic variables. The component S z is converted into 
the fermion occupation number and are creation and annihilation operators: 



S z n = *tttn-~, (3) 
Sn = exp I ITT ( S m + \) ) 



2 

(4) 



For the pure XXZ chain when A = 0, the ground state is a half-filled band with Fermi 
momenta kp = ±ir/2a where a is the lattice spacing. Low energy excitations around the 
Fermi surface can be described in terms of bosonic fields using a standard procedure AM : 



" " L 



(5) 
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Here L is the length of the system. The bosonic fields <&l,r describe left and right moving 
excitations. 

The random perturbation @, can also be expressed in terms of the bosonic fields and, 
at low energies, the continuum limit can be taken which results in an action pi H: 



-ii> 



S = J dxdr Qk(c\$) 2 + ^k(c\.$) 2 + g cos(2$) + rj{x)d x ^ + p(x)e^ + p*(x)e 

(6) 

for the field ^(x) = &r + Besides, g is a function of the anisotropy A, and r](x) 
and p(x) are composed of the Fourier components of the random field near zero and 2kp, 
respectively. Also, k is the spin stiffness of the ordered system given by the Bethe-ansatz 
solution Jl2[ ]: 

"-sC 1 -? -" 1 *)- (7) 

Afterwards, the replica formalism can be used to compute disorder averages of ground-state 
expectation values giving the effective action: 



s e 



77 = I dxdr 

a=l 

-J^D f dx dT dT 1 cos[$ a (x,r) - S^sc, r')], (8) 

a,/3 J 

where D =< (/if) 2 >= ^q/3. Thus, the disorder induces interactions between pairs of 
replicas. The interaction comes from averaging over p, that is the term due to the staggered 
component of the magnetic field. The term due to rj is irrelevant at large distances and is 
not included in (||). 

Using the fact that the dimension of the operator cos( < I ) ) is 1/4ttk, the RG eigenvalue 
A, which determines the rescaling of D is given by Q: 

A = 1 + 2z - 2C, (9) 

where the dynamical exponent z = 1 and C = l/(47r«). 

The sign of A determines the relevance of the randomness D. In the interval — 1/2 < 
A < 1, A is positive, indicating that the disorder is relevant and that the quasi-LRO is 
unstable under any small amount of randomness. On the other hand, A is negative in the 
region —1 < A < —1/2, so there the disorder is irrelevant and the quasi-LRO is expected 
to persist for small values of D. Therefore, in the latter case, a quantum phase transition 
between a quasi-LRO phase and a localized disordered phase has to occur for finite D Q. 
By means of a RG treatment for weak disorder Q the RG equations governing the critical 
properties can be obtained 0]. The flow lines in the k — D plane, near the point k = 1/67T 
are parabolas, (k — l/67r) 2 — D = being the critical one. The phase transition is of the 
Kosterlitz-Thouless type 0] and the correlation length corresponding to the spin correlation 
functions diverges like: 

£ ~ exp(,4/VA - A c ), (10) 

when the quasi-LRO phase is approached along a path in the phase diagram at fixed disor- 
der. In the region where the pure system is an Ising like ferromagnet (antiferromagnet) the 
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Figure 1: Phase diagram corresponding to the XXZ chain in a random field in the z direction 
according to . 



Imry-Ma argument [ 13 1 can be used to conclude that an arbitrarily weak random z field 
will destroy the long-range order. 

The phase diagram ho — A proposed by Doty and Fisher is shown in Fig. |l|. Note that 
the top of the region where the quasi-LRO persists i.e the maximum critical field is not 
established since the perturbative renormalization group result is valid only near A = —0.5. 

A numerical work analyzing this phase diagram was performed by Runge and Zimanyi 
||. They use an exact diagonalization procedure up to 16 sites and compute, among other 
things, the superfluid density p s which measures the sensitivity of the system to changes in 
the boundary conditions: 

1 d 2 E 



Ps = 7 a/nm2 ' ( U ) 



Here, 9 is a twist introduced at the boundary by adding to the Hamiltonian (|l|) the term 

' Si SI e ie + S^Sie-* 8 ) + A SfS z L . (12) 



Doing this, the chain is closed in such a way that a boson acquires a phase e w (e~ ie ) if it 
is hopping to the right (left) through the boundary. 

This superfluid density is the order parameter of the transition (p s > in the quasi-LRO 
phase and p s = in the localized phase). 

In H , the authors find a qualitative agreement with the RG predictions [||, [|, [| . However 
they conclude that it is necessary to study larger chains to obtain more accurate results 
near the critical region. 

In this paper we compute the superfluid density and the correlation length using the 
DMRG approach which allows us to study longer chains. In that way we expect to obtain 
further evidence and complementary results. 
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3 Numerical method 



We compute the properties of the system using an infinite-size DMRG algorithm 0. The 
disorder is included by doing the calculations for a single configuration of the random field 
and then, averaging over many realizations. Since the ground-state is non-magnetic in the 
quasi-LRO regime for the pure system and is expected to be also non-magnetic in the whole 
range of A when the disorder is turned on, we perform our calculations constrained to the 
Hilbert subspace where the z-component of the total spin is zero. We have checked that 
this approximation is justified after sufficient averaging is performed. 
We calculate the spin correlations 



C j (r) = (SZS? + S y Sy), (13) 
using open boundary conditions and the same superblock configuration B • »B as in |l 



(where B represents a block and • a single site). In (13), (. . .) denotes quantum expectation 
value in the ground state and j is a given realization of the random field. To avoid the 
boundary effects, we first grow the system until the length of each block is around 30 sites. 
After that we start to compute the correlations taking the spin operator Sq corresponding 
to the left site added at that DMRG step. In the next iteration this operator is rotated 
and truncated using the eigenvectors of the density matrix and so it will correspond to the 
right-most site of the left block. Iterating this procedure for the same operator we obtain 
after r DMRG steps an operator associated to a site which is at a distance r from the 
center. Simultaneously, the correlation Cj(r) is computed at each of those DMRG steps as 
(Sq + SqS~). Using this standard procedure it is only necessary to store one spin matrix 
(S^~). Afterwards, we average the logarithm of the correlation functions Cj(r) over disorder 
to obtain what is called the typical correlation [15|. Usually, more than 150 configurations 
of the random field are required. 

At large distances Cj(r) should be exponential outside the critical region where the dis- 
order is relevant (finite correlation length) and power law inside the critical region, where 
it is irrelevant (infinite correlation length). However, due to the truncation of the DMRG 
and the finite size effects, we also find an exponential behavior inside the critical region. 
This fact makes possible to define a typical correlation length from the slopes of the curves 
In Cj(r) vs r and analyze its behavior in the phase diagram A — ho. The correlation length 
computed in this way is of course approximated and consequently gives qualitative infor- 
mation about the phase diagram. 

If chains of various length are considered and, if for each of them, various number of 
states are kept, an extrapolation to infinite length and infinite number of DMRG states 
can be performed to decrease the error. This is discussed in section [| where we do the 
extrapolations for a few special points in the phase diagram. 

To compute the superfluid density which is the order parameter of the transition between 
the quasi-LRO and the localized states, we need to consider a closed system with twisted 
boundary conditions. Closed systems can also be treated with the density matrix algorithm 
as shown in [ 16 1 . In that work two superblock configurations were analyzed: B • B* and 
B • »B. It was found that the first one is better since it has a smaller truncation error 
(because the two blocks are not connected directly), although larger than the error obtained 
with open boundary conditions. In our case, we found the configuration B • »B more 
convenient to introduce twisted boundary conditions since the twist can be introduced 
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Figure 2: Superblock configuration used for the calculation of the superfluid density. The 
twist is introduced between the two single sites added at each step of the DMRG. 



between the two single sites that we add at each step of the procedure (in that case the 
operators S + and S~ which appear in ( |i~2| ) are not truncated). On the other hand since we 
are also interested in minimizing the truncation error, for the actual calculation we chose 
the configuration in Fig. [2] which combines the advantages of both. Note that the left and 
right blocks have different number of sites although the same number of states are kept on 
both of them. This asymmetry implies that we have to compute two density matrices, one 
for the left block and one for the right one. However this is necessary anyway since the 
reflection symmetry is already broken by the random magnetic field. 

To approximate eq. ([ll]) numerically, given a configuration of the random magnetic 
field, we perform the DMRG procedure twice (with and without twist) and compute for 
each length L the finite difference 



Ps{L) 



1 E(LV9) - E(0) 



L 



(V9f 



(14) 



where we used that the first derivative dE/dV9 vanishes at V# = due to the symmetry 
E(X78) = E(— V#). Afterwards, we average this quantity over many configurations of the 
field. The h larger the strength of the field the more configurations needed (in most of the 
cases more than 200 configurations were considered). 

To test the convergence of the DMRG method we computed p s keeping different number 
of eigenstates of the density matrix. It was found that the DMRG procedure converges if 
we keep V# constant i.e. the twist 6 increases with the length of the chain. This is shown 
in Fig. ||, where p s vs. L for A = —0.75 is plotted for a random field ho = 0.1 keeping up to 
128 eigenstates of the density matrix. Since we need LS/9 to be small, in the calculations 
we took W = 0.01 and grew the system up to 50 sites. 

We also checked (see Fig. that without disorder our numerical results for p s agree 
with the exact solution 

7T sin p 



Ps (A,h = 0) 



2p(ir - pY 



(15) 



found in fi~0[|. Here, p 



cos 
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Figure 3: The dependence of the superfluid density with the length of the chain for A = 
—0.75 and a random field with strength ho = 0.1 keeping 16, 32, 64 and 128 eigenstates of 
the density matrix. It is possible to observe the convergence of the method. 
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Figure 4: The dependence of the superfluid density with the anisotropy for the system 
without disorder. Our numerical results for a chain of 50 sites keeping 32 eigenstates of the 
density matrix are compared to the exact solution found in [fil]] . 
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Finally, since we are truncating the Hilbert space it is important to estimate how good 
this approximation is. A standard measure of that is the truncation error defined as the sum 
of the eigenvalues of the density matrix that are discarded. Keeping up to 32 eigenstates 
of the density matrix we obtained a truncation error less than 10~ 7 for open boundary 
conditions (spin correlations) and less than 10 -5 for twisted boundary conditions (superfluid 
density) . 

4 Results 

4.1 (S?S? + SfS^) correlation 

Before summarizing the correlation functions results, let us discuss what should be the 
behavior of the correlation length in a phase diagram like the one proposed in Fig. |l[ As 
we said before, in the localized phase it is possible to define a correlation length, i.e. the 
correlations functions display an exponential behavior at large distances. Instead, in the 
phase with quasi-LRO this is not possible because, since the region is critical, the correlation 
length is infinite. This means that, if we were able to solve the problem exactly, we should 
observe in a 3-dimensional plot of correlation length vs (A,/jq) that this quantity diverges 
when we approach the line that divides the two regimes. This is not the case when we 
consider a finite chain and moreover when we truncate the Hilbert space. What we should 
obtain instead is that the correlation length has a peak located in a region that coincides 
approximately with the critical region. This peak should become sharper, higher and closer 
to the transition as we consider larger chains and keep more DMRG states. 
Having this in mind, we computed the spin-spin correlations 

0(L) = In | (S§Si + S%Sl) I, (16) 

transverse to the random field for several points of the phase diagram A — /to- Here, the 
bar indicates average over disorder. 

In Fig. |5|, Ct(L) is plotted for ho = 0.3 and different values of A. The linear behavior 
at large distances means the correlations decay exponentially. 

From the slope of these lines we estimate the typical correlation length £ and analyzing 
its behavior as a function of A and ho we obtain the 3-dimensional plot shown in Fig. |6| for 
16 and 32 DMRG states. There, we can see that the correlation length not only increases 
when we approach the axis ho = (pure system) keeping A constant but also in a region 
with finite disorder centralized in a negative value of A between —1 and —0.5. This is an 
indication of the presence of a critical region for finite disorder where the quasi-LRO found 
for ho = persists. 

To obtain an approximation to the shape of the critical region, we propose to analyze the 
contour plot of the typical correlation length. Following the previous discussion, in the limit 
of infinite length and infinite number of DMRG states, the contour lines will concentrate 
on the border of the critical region. If, for a given truncation, we use the input that the 
critical region starts at A = —0.5, the choice of a contour line such that ho ~ at A = —0.5 
will give us qualitative information about the shape of the critical region. This choice will 
be justified in the next section when we discuss the superfluid density results. Fig. ^ shows 
the contour plot of the typical correlation length for 16 and 32 DMRG states. From the 
chosen contour lines (bold lines in the figure) we estimate /t C rit.(A = —0.75) ~ 0.17 for 16 
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Figure 5: Average of the logarithm of the correlations (typical correlations) vs. the distance 
for ho = 0.3 and different values of the anisotropy. The statistical error originated in the 
random average is, in all the cases, smaller than the size of the symbols. 

states kept and /i cr ;t.(A = —0.75) ~ 0.20 keeping 32 states. We want to remark here that 
in Fig. 0, the contour lines used to estimate the boundary of the critical region correspond 
to values of £ smaller than 100 (see Fig ^) while the chains considered were larger than 200 
sites. Thus, the finite size effects are small and the main approximation is in the truncation 
of the Hilbert space. 



4.2 Superfiuid density 



We computed the superfiuid density using (14) and performed afterwards the average over 
disorder. In Fig. ||[ we show p s vs L for different values of the strength of the random 
field ho and A = —0.75. As expected, for a fixed length of the chain, p s decreases with ho. 
Note that we obtain non zero values of p s because we study finite systems. For the infinite 
system the superfiuid density is different from zero only for (ho = 0,|A[ < 1) and inside 
the critical region, that is, where the correlation length diverges. For that reason p s can be 
used to study the phase diagram of Fig. |l]. 

First we are going to check that the quasi-LRO region extends up to A = —0.5. This 
point is characterized by the fact that there, the correlation length critical exponent diverges. 
In the region —1/2 < A < 1 where the disorder is relevant, the dependence of the correlation 
length with the disorder D = /iq/3 near the critical line corresponding to ho = 0, should be 
£ ~ D~^ a and according to the perturbative renormalization group results [||, the critical 
exponent <j) s diverges at A = —0.5 like 

(17) 



3 — (1 — — arccos A) 1 

To extract this exponent from our numerical data, we analyzed the derivative dp s /dD at 
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Figure 6: The dependence of the typical correlation length in the plane A — ho. a) Keeping 
16 eigenstates of the density matrix, b) Keeping 32 eigenstates of the density matrix. 
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Figure 8: The dependence of the superfluid density with the number of sites for different 
values of the random field and A = —0.75. The statistical error originated in the random 
average is, in all the cases, smaller than the size of the symbols. 



D = 0. Assuming that p s = p s (L/£) for L ^ 1 ||, the behavior of ln(dp s /dD\D=o) with 
InL should be linear with the slope corresponding to <f)~ . In the inset of Fig. ^ we show 
this linear behavior for different values of A. The critical exponent extracted from our data 
and the analytical prediction (eq. |l7|) are shown in Fig. ^ where a good agreement can be 
observed. This result is in line with the existence of a critical region for A < —0.5 where 
the disorder is irrelevant and justify our choice made in the £ contour plots. 

For a given truncation, the border of this region can be estimated studying the de- 
pendence of p s with the strength of the random field for different lengths of the chain, L. 
Figure 1C is a plot of p s vs ho for A = —0.75 and different lengths. The exact solution 
(infinite length and infinite DMRG states) should show a discontinuity in a transition of 
the Kosterlitz-Thouless type {p s {ho < /i cr it.) > and p s (ho > fa C rit.) = 0). 

Using this fact, an estimation of the critical field can be obtained analyzing the behavior 
of the function 5p s (A, ho) = p s (A,ho, L) — p s (A,ho,L') where L and L' are fixed and 
L,L' 3> 1. This corresponds, in Fig. to taking two of the curves and computing the 
difference. If the two curves correspond to lengths that are L, L' S> 1, the difference between 
them will start to increase considerably only when ho > h CT it_. This criterion assumes that 
for L 1 p s (L,£) = p s (L/£) and then, if the truncated system correlation length is much 
larger than L, p s becomes almost independent on the length L. In this way, this estimation 
is affected mainly by the truncation of the DMRG method and then, we expect to have 
8p s almost constant in the region where the truncated system correlation length is much 
larger than L = 50 (which is the maximum length considered in the supernuid density 
calculations). In Fig. |ll] we show 5p s (A,ho) vs. ho (taking L = 30 and V = 50) for 
A = —0.5 and A = —0.75. The points in these curves are obtained by averaging over 
disorder and have a statistical error which is shown in the figure. As it is possible to 
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Figure 9: The dependence of the correlation length critical exponent with the anisotropy. 
The stars correspond to our numerical data and the solid line is the perturbative renormal- 
ization group prediction [II. The inset shows the dependence of ]n(dp s /dD\i)=o) with InL 
for different anisotropy values; the slopes correspond to the inverse of the correlation length 
critical exponent. 
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Figure 10: The dependence of the superfluid density with the strength of the random field 
for systems with different number of sites. The the anisotropy value is A = —0.75. 
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Figure 11: The dependence of 5p s with ho for A = —0.5 and A = —0.75. The error bars 
indicate the statistical error originated in the random average. 



observe, for A = —0.75 the function stays almost constant until a value of ho ~ 0.25. We 
have checked that this behavior is the same taking different values of L and V . 

We also analyzed 5p s (A, ho) for other anisotropy values. The results are shown in Fig. 



12| . Remarkably, all the curves corresponding to values of A between —1 and —0.5 fall 
between those corresponding to A = —0.5 and A = —0.75 supporting the statement that 
the maximum critical field occurs for A = —0.75. An estimation of the critical field is then 
given by the intersection of the curve with the /io-axis. To obtain that intersection we use 
the first two points which differ from zero by more than the statistical error. Also from the 
error in these points we estimate the error in the critical field so computed. 



In Fig. 13, we plot these values of the critical field in the phase diagram giving an 



approximation to the boundary of the critical region. 



5 Extrapolations for A = —0.75 

In this section we study the special case of A = —0.75 and perform the extrapolations to 
infinite length and infinite number of DMRG states of the computed properties for selected 
values of the random magnetic field. 

For the superfluid density p s , we considered two points on the phase diagram: (A = 
—0.75, ho = 0.1) and (A = —0.75, ho = 0.5). From the approximated phase diagram of 
Fig. ^ we expect ho = 0.1 to be inside the critical region and ho = 0.5 to be outside. 



In Fig. |b| we plot p s as a function of 1/L. At large L the two cases display opposite 
behavior. For ho = 0.1, p s increases as we keep more DMRG states whereas for ho = 0.5, 
it decreases. This is already an indication that the point (A = —0.75, ho = 0.1) is inside 
the quasi-LRO region whereas the point (A = —0.75, ho = 0.5) is outside. To find the 
asymptotic values po(s) = p s (s, L — > oo), we fitted each curve in Fig. [14| using the expression 
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Figure 12: <5/5 s vs /iq for different values of A keeping 32 eigenstates of the density matrix. 
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Figure 13: Values of the critical field (keeping 32 states) obtained from Fig. 12 as explained 
in the text. The dashed curve corresponds to the phase boundary estimated by the analysis 
of the correlations. 
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p s (s,L) = po(s) + a(s) / L@( s \ The solid lines in Fig. 14 are the results of these fittings. 
Finally, the values of po(s) vs. 1/s are plotted in Fig. 15. Again, we considered the 
expression po(s) = p^ + a/s b to fit the curves. For ho = 0.1 the extrapolated value is 
Poo = 0.586 ± 0.003, confirming that the point (A = —0.75, ho = 0.1) is inside the critical 
region whereas for ho = 0.5 the best fit was obtained for poo = 0, indicating the point 
(A = —0.75, ho = 0.5) is in the localized phase as we expected from the previous analysis. 



We proceeded in the same way with the correlation length. In Fig. 16 we plot the 
typical correlation length £ as a function of 1/L for A = —0.75 and ho = 0.1, ho = 0.3 and 
ho = 0.5. In the ho = 0.1 case, the correlation length increases drastically as we keep more 
and more DMRG states. As more states are kept, the larger has to be the chain to reach the 
asymptotic value corresponding to the infinite system. Instead, for ho = 0.3 and ho = 0.5, 
the behavior is different. Not only the asymptotic values for the correlation length are one 
order of magnitude smaller but also the curves become closer as we increase the number of 
DMRG states indicating that the correlation length is finite. As we did for the superfluid 



density, we found the asymptotic values £o(s) = £( s > L — > oo), fitting each curve in Fig. 16 
with the expression f(a,L) = f„(a) + a{s)/L^ s \ The lines in Fig. [fij are the results of 
these fittings. The asymptotic values £a(s) = £(s,L — > oo) found in this way are plotted as 
a function of 1/s in Fig. [H]. In all the cases, the errors in the data are smaller than 1% 
and then, smaller than the size of the symbols. 



To fit the curves in Fig. [17| we used again the expression £o(s) = £oo + a/s . For 



ho = 0.1 the best fit was obtained for a negative exponent (b = —0.155) implying an infinite 
correlation length for infinite number of DMRG states. On the other hand, for ho = 0.3 
and ho = 0.5 we find b possitive and the extrapolated values for the correlation length are 
^(ho = 0.3) = 108 ± 18, and ^(ho = 0.5) = 26.45 ± 2.95. 

In summary, from the extrapolations performed for A = —0.75, we conclude the exis- 
tence of a quasi-LRO and put bounds on the exact value of the critical field. This value 
should be larger than ho = 0.1 and smaller than ho = 0.3. 



6 Conclusions 

In this work we analyzed numerically the quantum phase diagram of the spin 1/2 XXZ 
chain in a random magnetic field pointing in the z direction. 

The two-point spin-spin correlation functions in the XY plane and the superfluid den- 
sity were computed using an infinite-size density matrix renormalization group (DMRG) 
algorithm. 

The asymptotic exponential behavior of the spin correlations allowed us to extract the 
typical correlation length and study its dependence with the anisotropy A and the strength 
of the random magnetic field ho- The contour plots of the correlation length as a function 
of A and ho give numerical evidence of the existence of a critical region for finite disorder 
and negative values of the anisotropy. In this region the correlation length increases as we 
keep more states in the DMRG method. Using our data for p s we extracted the behavior of 
the correlation length critical exponent. A good agreement with the analytical prediction 
was found Q supporting the idea that the critical region starts in A = —0.5. Then, an 
approximation to the boundary of the critical region can be obtained from the contour plots 
of the typical correlation length taking a contour line such that for A = —0.5 the value of 
the field is ho — 0. In this way, we estimated /i cr i t .(A = —0.75) ~ 0.20 for 32 states kept. 
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Figure 14: The Superfluid Density for A = —0.75 as a function of the inverse of the number 
of sites for different number of states kept in the DMRG. a) ho = 0.1. b) ho = 0.5. The 
solid lines are the fittings obtained using the expression p s (s,L) = p (s) + a{s)/L^ s \ The 
statistical error originated in the random average is smaller than the size of the symbols. 
At large L the two cases display opposite behavior as we keep more DMRG states. 
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Figure 15: Extrapolated values po as a function of the inverse of the number of DMRG 
states. The lines are the result of the fittings using the expression po(s) = Poo + a/s b . For 
ho = 0.1 the extrapolated value is p^ = 0.586 ±0.003, whereas for ho = 0.5 the best fit was 
obtained for p^ = 0. The error bars were estimated from the statistical errors of the data. 
In the ho = 0.1 case they are smaller than the size of the symbols. 

From the analysis of the function 5p s (A, ho) = p s (A, ho, L) — p s (A, ho,L') for L, V 3> 1, we 
obtained an independent estimation of the critical field /i C rit. (A) for several values of A and 
drew the corresponding points of the line that separates the quasi-LRO regime from the 
localized one. We found that the largest value of the critical field occurs at A = —0.75 and 
is approximately /i cr i t .(A = —0.75) = 0.25 for 32 states kept, in close agreement with the 
results obtained from the correlation length analysis. Since these results are approximated 
and are mainly affected by the truncation of the Hilbert space, to further support them, 
we extrapolated the typical correlation length and the superfluid density to infinite length 
and infinite number of DMRG states for the special case of anisotropy, A = —0.75, and 
selected values of the random magnetic field. The extrapolations allow us to conclude the 
existence of a quasi-long-range order region and put bounds on the critical field. The range 
of possible values for the critical field (0.1 < h cv \t. < 0.3) contains the approximated results 
before obtained. 
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Figure 16: The typical correlation length for A = —0.75 as a function of the inverse of 
the number of sites for different number of states kept in the DMRG. a) ho = 0.1. b) 
ho = 0.3 c) ho = 0.5. The lines are the fittings obtained using the expression £(s,L) = 
£o(s) + a(s)/L /3 ^ s \ The statistical error originated in the random average is smaller than 
the size of the symbols. 
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Figure 17: Extrapolated values £o as a function of the inverse of the number of DMRG 
states. The lines are the result of the fittings using the expression £o( s ) = Coo + a/s b . 
The error bars are smaller than the size of the symbols. For ho = 0.1 we find b negative 
implying an infinite correlation length for infinite number of DMRG states. For ho = 0.3 
and ho = 0.5 we find b possitive and the extrapolated values for the correlation length are 
^(ho = 0.3) = 108 ± 18 and ^(ho = 0.5) = 26.45 ± 2.95. 
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